Multidimensional solitons in a low-dimensional periodic potential 
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Using the variational approximation(VA) and direct simulations, we find stable 2D and 3D 
solitons in the self-attractive Gross-Pitaevskii equation (GPE) with a potential which is uniform 
in one direction (z) and periodic in the others (but the quasi- ID potentials cannot stabilize 3D 
solitons). The family of solitons includes single- and multi-peaked ones. The results apply to Bose- 
Einstein condensates (BECs) in optical lattices (OLs), and to spatial or spatiotemporal solitons in 
layered optical media. This is the first prediction of mobile 2D and 3D solitons in BECs, as they 
keep mobility along z. Head-on collisions of in-phase solitons lead to their fusion into a collapsing 
pulse. Solitons colliding in adjacent OL-induced channels may form a bound state (BS), which then 
relaxes to a stable asymmetric form. An initially unstable soliton splits into a three-soliton BS. 
Localized states in the self-repulsive GPE with the low-dimensional OL are found too. 
PACS numbers: 03.75.Fi, 05.30.Jp, 05.45.-a 



Introduction. Solitons in multidimensional nonlinear 
Schrodinger equations (NLSEs) with a periodic potential 
have recently attracted considerable interest. In particu- 
lar, self-trapping of spatial beams in nonlinear photonic 
crystals is described by a 2D equation. In this case, de- 
spite the possibility of collapse [1], simulations reveal ro- 
bust 2D solitons in the self- focusing model [2] . A similar 
medium can be created by a grid of laser beams illumi- 
nating a photorefractive sample [3]. 

Similar 2D and 3D models with a periodic potential 
describe a Bose-Einstein condensate (BEC) trapped in 
an optical lattice (OL) [4]; this application is relevant, 
as experimental techniques for loading BECs into mul- 
tidimensional OLs were recently developed [5]. Stable 
solitons can be supported by an OL even in self-repulsive 
BECs [6,7]. In the case of self-attraction, 2D and 3D soli- 
tons (including 2D vortices) are stable in the self-focusing 
model with the OL potential [8], despite the possibility 
of the collapse. 

An issue of direct physical relevance, and a subject of 
the present work, are multidimensional solitons in media 
equipped with periodic potentials of a lower dimension, 
viz., quasi-lD (Q1D) and Q2D lattices in the 2D and 3D 
cases, respectively. In optics, the 2D equation in the spa- 
tial domain governs the beam propagation in a layered 
bulk medium along the layers, which is an extension of 
a ID multichannel system introduced in Ref. [9], where 
the potential was induced by transverse modulation of 
the refractive index (RI). In the temporal domain, the 
2D and 3D equations govern, respectively, the longitu- 
dinal propagation of spatiotemporal optical solitons in a 
layered planar waveguide, or in a bulk medium with the 
RI periodically modulated in both transverse directions. 
The 2D and 3D cases directly apply to BECs loaded in 
a Q1D or Q2D lattice. The physical significance of this 
setting is two-fold: first, in the experiment it is much 



easier to create a lower-dimensional lattice than a full- 
dimensional one, both in BECs and in optics, hence the 
use of such lattices is the most straightforward way to 
create multidimensional solitons in these media. Second, 
the solitons created this way can freely move in the un- 
confined direction, which also suggests a possibility to 
study their collisions, and to look for their bound states 
(BSs). As yet, no other way to create multidimensional 
mobile solitons in BECs and their BSs has been proposed. 
In optics, solitons of this type suggest new applications. 
Indeed, in an optical medium with the full-dimensional 
periodic potential, transfer of a trapped beam from one 
position to another is difficult, as the necessary exter- 
nal push strongly disturbs the beam [9]. In the lower- 
dimensional potential, the beam can slide along the free 
direction, making the transfer easy. In BECs confined 
by a low-dimensional OL, matter- wave solitons can be 
driven in the free direction by a weak laser beam. 

The model and VA. The normalized form of the self- 
focusing 2D NLSE with a Q1D periodic potential of the 
strength e is 



iut + u xx + u yy + [ecos(2x) + xM' 



0. 



(1) 



where x = ±1 in the case of the self-attraction/repulsion. 
In BECs, t is time, while in optics it is the propagation 
distance. For BECs or spatial optical solitons, x and 
y are transverse coordinates; for spatiotemporal optical 
solitons in a 2D waveguide with anomalous chromatic 
dispersion, y is a temporal variable. The 3D version of 
Eq. (1) is 

iut + u xx + Uyy + u zz + {e [cos(2a;) + cos(2y)] + x\u\ 2 } u = 0. 

(2) 

Equations (1) and (2) conserve the Hamiltonian, the 
norm JV = J|u(r)| dr (the number of atoms in the BEC, 
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or total power/energy of the spatial/spatiotemporal opti- 
cal soliton), and the momentum along the free direction. 
The equations are normalized so that the period of the 
potential is n, the free parameters being e and N. 

Stationary solutions to Eq. (1) are looked for as 
u(x,y,t) = U(x, y) exp(— i/J,t), with a chemical potential 
/j, (propagation constant, in optics), which yields 

fiU + U xx + U yy + [e cos(2a;) + X U 2 ] U = 0, (3) 

and similar in the 3D case. To apply the VA 
(variational approximation) to Eq. (3), we adopt 
the ansatz U = A exp [— (as 2 + by 2 ) /2] , or U = 
A exp [- (a(x 2 + y 2 ) + bz 2 ) /2] in the 3D case, with the 

norms N 2 u = nA 2 /\/ab, and A 3D = n 3 / 2 A 2 /(aVb). Fol- 
lowing the procedure of the application of VA to multi- 
dimensional solitons [10], we derive variational equations 
from the Lagrangian of Eq. (3). In the 2D case, they are 
(N ee er x l a ) 

N X = {At: /a) V 'a 2 - 2eH, (j, = -a - eN(l - 3/a), (4) 
and in the 3D case [in both cases, b = (a 2 — 2eH) / a] , 

N X = 2 (27r/a) 3/2 V 'a 2 - 2ett, (j, = -a/2 - eK (2 - 3/a) . 

(5) 

Evidently, solutions are possible only if \ > (self- 
attraction). The 2D solutions exist in the interval 

A min ee 4vr v'l - e/e C r < N < V max ee 4tt (6) 

[see Fig. 1(a)], if e < e cr ee e 2 /8 w 0.924 (A min = if 
£ > £ cr ), while 3D solutions can be found for any N. The 
fixed lattice period rules out establishing general scaling 
relations for these soliton families; in the case of very 
large e, the potential valleys become isolated, then one 
is actually dealing with cigar-shaped [11] Q1D solitons, 
that obey the usual ID NLSE scaling. 

The solution families which definitely meet the condi- 
tion dfi/dN < are expected to be stable pursuant to 
the Vakhitov-Kolokolov (VK) criterion [12], while ones 
with d/i/dN positive or nearly zero should be unstable. 
It follows from Eqs. (4) and (5) and is evident in Figs. 
l(b,c) that the 2D solitons are stable in the existence in- 
terval (6). At small e, the 3D solutions arc unstable for 
any N. Equation (5) predicts that a stability interval of 
a width A./V ~ yfe — £q appears around jV « IO.Itt 3 / 2 if 
e exceeds eo « 0.242. 
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FIG. 1. The numerically found (joined squares) and 
VA-predicted (dashed lines) existence limits for stable 2D soli- 
tons (a). VA-predicted dependences y,(N) for the 2D (b) and 
3D (c) solitons, which determine their VK stability. 

Numerical results. Solitons are generated by the 
imaginary-time evolution method [13], starting with the 
VA-predicted waveforms. Stability is verified by direct 
simulations of perturbed solitons in real time. Typi- 
cal examples of the thus found stable 2D and 3D soli- 
tons are displayed in Fig. 2 and 3, the 3D ones be- 
ing nearly isotropic in the (x, y) plane and elongated 
in the free direction z. As is seen in Fig. 1(a), the 
existence limits for the stable 2D solitons are close to 
the VA prediction, unless N is small; in that case, the 
discrepancy is due to the fact that the soliton develops 
a multi-peaked structure [Fig. 2(b)], which is not ac- 
counted for by the simple ansatz adopted above. The 
ansatz can be generalized for this case, setting U — 
(A + B cosx) exp [— (as 2 + by 2 ) /2] , but the final result 
is then quite messy. A numerically found stability re- 
gion for the 3D solitons is also quite close to the VA 
prediction, and multi-peaked 3D solitons relate to their 
single-peaked counterparts (Fig. 3) the same way as Fig. 
2(b) to Fig. 2(a). 




FIG. 2. Examples of single- and multi-peaked stable 2D 
solitons in the Q1D potential (e = 2). The norm of the soli- 
ton is N = 2n (a) and N = (b). 

The VA predicts that, in the 3D case with the Q1D 
(rather than Q2D) lattice, all the solitons are VK- 
unstable. Accordingly, simulations never produced stable 
solitons in this case. This feature can be explained by the 
fact that, in the free 2D subspace, the soliton is essen- 
tially tantamount to the unstable Townes soliton [1]. In 
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the 2D NLSE with the 2D lattice, stable solitons with in- 
trinsic vorticity were also found in Ref. [8] . Vortices were 
found in the present 2D model too, but they are always 
unstable. 




FIG. 3. A stable single-peaked 3D soliton in the Q2D po- 
tential (N — 2tt, e — 5.0) is shown through its z = (a) and 
y — (b) cross sections. 

We have also investigated the case of self-repulsion, 
X = — 1, when the low-dimensional lattice potential 
cannot support a completely localized pulse. However, 
adding a usual parabolic trap readily gives rise to stable 
states, which are solitons across the lattice and Thomas- 
Fermi states along the free direction. It is noteworthy 
that, with x < 0, these states assume a multi-peaked 
or single-peaked shape under the action of a weaker or 
stronger OL potential, respectively, which is opposite to 
what was reported above for x > 0, cf. Fig. 2. An 
explanation is that, in the case of the self- repulsion and 
e = 0, no solitons exist (even unstable ones, like the 
above-mentioned Townes soliton). 

The influence of the parabolic trap was checked too in 
the case of \ = +1. If the 2D or 3D soliton is displaced 
from the central position, it performs harmonic oscilla- 
tions along the free direction, completely preserving its 
integrity. Actually, the mobility in the free direction is 
the most essential difference of the present multidimen- 
sional solitons from those predicted in other BEC models 
[6-8]. This suggests new possibilities, such as collisions 
between moving solitons. 

Simulations demonstrate that 2D and 3D in-phase soli- 
tons which collide head-on with velocities ±v, which are 
below a critical value v c , merge into a single pulse, whose 
norm exceeds the critical value iV max [see Eq. (6)], hence 
it quickly collapses. If v > v c , the solitons pass through 
each other (for instance, v c = 8.5 for e = 2.0, when 



the norm of each 2D soliton is N — 2.bn), which is ex- 
plained by the fact that the collision time, ~ is then 
smaller than the collapse time, ~ 1/iV. Colliding 7r-out- 
of-phase solitons always bounce back, and two such soli- 
tons, placed inside the trap, perform stable oscillations 
with periodic elastic collisions. 



In the low-dimensional potential, collisions are also 
possible between solitons moving in adjacent "tracks" 
(channels). In that case, they pass each other quasi- 
clastically if the collision velocity is large. If it is not 
too large, each soliton captures a small "satellite" in the 
other channel. If the velocity falls below a critical value 
for given e, the outcome of the collision is altogether dif- 
ferent: two solitons form a quiescent symmetric bound 
state (BS) [Fig. 4(a)]. Later, the symmetric BS develops 
an intrinsic instability, spontaneously rearranging itself 
into a completely stable asymmetric BS [Fig. 4(b)]. This 
behavior resembles the symmetry-breaking instability of 
optical solitons in dual-core fibers [14] and the macro- 
scopic quantum-self-trapping effect in BECs confined to 
a double- well potential [15]. In fact, this is the first exam- 
ple of a stable BS of two solitons in BECs. Formation of 
more complex stable 2D and 3D BSs also occurs in a dif- 
ferent way: VK-unstable waveforms predicted by the VA 
undergo violent evolution, shedding ~ 50% of their norm 
and eventually forming a BS of three solitons, with a nar- 
row tall one in the middle, and small-amplitude broad 
satellites in the adjacent channels [sec Fig. 5]. 




FIG. 4. A symmetric bound state (BS) of 2D in-phase soli- 
tons formed in adjacent tracks as a result of collision between 
them (a), and its subsequent spontaneous rearrangement into 
a stable asymmetric BS (b). The initial norm and velocities 
of the solitons are N — 2tv and v = ±0.1, and e = 3. In the 
3D case, the scenario is quite similar. 
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FIG. 5. Bound state of three solitons developed 
from the VK - unstable 2D soliton with parameters 
A = 0.63, a = b = 0.1, and e = 10. Half of the initial 
norm is lost during the formation. 

Conclusion. We have demonstrated that the peri- 
odic potential whose dimension is smaller by 1 than the 
full spatial dimension readily supports stable single- and 
multi-peaked solitons in 2D and 3D self-focusing media 
(although the quasi-lD potential cannot stabilize 3D soli- 
tons), which suggest new settings for experimental search 
of solitons. 2D solitons exist in a finite interval of the 
values of the norm, and 3D ones are stable if the lat- 
tice strength exceeds the minimum value. In the case 
of self-repulsion with a parabolic trap, stable localized 
states are found too. The dependence of their structure 
on the lattice strength is opposite to that in the case of 
the self-attraction: with the increase of the strength, a 
single-peaked state is changed by a multi-peaked one. 

These solitons are the first example of mobile multi- 
dimensional pulses in BECs. Head-on collisions between 
the in-phase solitons may lead to their fusion and col- 
lapse, while out-of-phase solitons collide elastically in- 
definitely many times. Collision between solitons in ad- 
jacent tracks may create a bound state (BS), which then 
relaxes to an asymmetric stable shape. Three-soliton BSs 
are created as a result of the evolution of initially unsta- 
ble solitons. 
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